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Abstract. Sample type and preservation methods are likely to influence the results of microbiome analyses. Here, we sam- 
pled tadpoles of Nanorana parkeri non-invasively for feces, and subsequently for hindgut through dissection. Gut and fe- 
cal samples were stored in ethanol, and additionally, part of the gut samples were frozen, but temporarily thawed during 
transport as it often happens under difficult field conditions. Our results showed that both substrate (gut vs. feces) as well 
as preservation method can influence the analysis of intestinal microbiomes. Frozen gut samples strongly differed from 
ethanol-stored samples, and most of the frozen samples were dominated (in relative abundance) by a set of Proteobacteria 
OTUs that were completely absent from the ethanol-stored samples. This blooming of contaminant bacteria occurred after 
less than 12 h of thawing, thus caution should be taken when constancy of cold temperatures cannot be maintained in the 
field for sample preservation purposes. Among ethanol-stored samples, whereas the recovered bacterial richness and phy- 
logenetic diversity did not differ, the community structure varied significantly. Thus, for meaningful gut microbiome me- 
ta-analyses it is important to acknowledge the effect of the sampling substrate for the microbial community composition. 
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Introduction 


Many studies have targeted the composition, dynamics and 
highlighted the health relevance of the human gut micro- 
biome (e.g. CLEMENTE et al. 2012, TUDDENHAM & SEARS 
2015). Likewise, animal-associated microbial communi- 
ties play important roles in the biology and health of their 
hosts (McFALL-Neat et al. 2013), including food degrada- 
tion and energy harvest (TREMAROLI & BACKHED 2012), 
immunity regulation (ROUND & MAZMANIAN 2009), and 
physical development (SOMMER & BACKHED 2013). These 
studies are typically based on high-throughput sequencing 
(HTS) of short 16S rRNA gene amplicons, where signifi- 
cantly different composition and diversity patterns of the 
host-associated microbiomes are driven by host taxonomy, 
ecology and environment (ENGEL & MORAN 2013, KOHL et 
al. 2014, SHAPIRA 2016, VENCES et al. 2016, LYRA et al. 2018, 
REESE & DUNN 2018). 

While the intestinal microbiome of mammals has been 
the subject of numerous studies (e.g. Li et al. 2016a, 2016b, 
SONNENBURG et al. 2016), research on microbiomes as- 
sociated with amphibians is still in its infancy. Work on 
this vertebrate group has typically focused on the cutane- 


ous microbiome (McKENZIE et al. 2012, KUENEMAN et al. 
2019), often in the context of its effects on pathogenic fungi 
causing amphibian declines (BLETZ et al. 2013, WALKE & 
BELDEN 2016). In particular, work on the amphibian gut 
microbiome is scarce (BLETZ et al. 2016, WENG et al. 2016, 
KNuTIE et al. 2017, FONTAINE et al. 2018, HUANG et al. 2018, 
Lyra et al. 2018, WIEBLER et al. 2018) but of substantial bio- 
logical interest especially in frogs, given the drastic restruc- 
turing of the gut during metamorphosis from a largely her- 
bivorous tadpole to an exclusively carnivorous frog (KOHL 
et al. 2013, VENCES et al. 2016). 

Understanding which factors influence microbiomes 
at local and global scale is of importance to unravel gen- 
eral biogeographic and macro-ecological trends (AMEND 
et al. 2013, BAHRAM et al. 2018b, KUENEMAN et al. 2019) 
and host-microbiome interactions (REESE & DUNN 2018), 
including diseases. Sampling the microbiome of the intes- 
tine is a moderately to strongly invasive approach, espe- 
cially in small animals where it requires killing the studied 
individuals. Instead, samples of fresh feces are often used 
and the fecal microbiome is considered as representative 
of the gut microbiome (e.g. REESE & DUNN 2018, ZIERER 
et al. 2018). However, indications exist that fecal character- 
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istics influence the composition and richness of detected 
microbiota (FALony et al. 2018). Therefore, it is uncertain 
to which extent the fecal microbiome may serve as proxy 
for the gut microbiome across hosts, and whether commu- 
nities from feces and gut samples are fully comparable (e.g. 
ERICSSON et al. 2018). In particular, only few studies ad- 
dress this question in tadpoles (Sonc et al. 2018). 

For large-scale analyses, especially meta-analyses of 
data sets originating from a diversity of sources, compa- 
rability of data is a basic prerequisite. This refers not only 
to laboratory methods but also extends to sampling and 
sample preservation (FRICKER et al. 2019). Thus, a further 
factor influencing the inference of microbiome composi- 
tion from high-throughput sequencing of amplicons is the 
method of sample preservation. Freezing samples at -20°C 
immediately upon collection has been defined as the gold 
standard to ensure the microbial community does not 
change until DNA extraction (SONG et al. 2016). Keeping 
samples uninterruptedly at this temperature is however 
not always possible under difficult field conditions, which 
might lead to alterations of the microbial community com- 
position during episodes of thawing (CHoo et al. 2015). 

In our study, we tested the effect of sampling substrate, 
feces vs. hindgut, to characterize the gut microbiota of tad- 
poles of the High Himalaya frog (Nanorana parkeri, Di- 
croglossidae). In addition, we compared the similarity of 
the detected microbiota using freezing vs. ethanol (EtOH) 
preservation method of the samples from the same spec- 
imens. The sample types per specimen included EtOH- 
stored feces, EtOH-stored gut and frozen gut samples. As 
feces samples are widely considered to reflect gut micro- 
biota, we expected similar richness and community pat- 
terns from the EtOH-stored feces and gut samples. In this 
study, the frozen gut samples were exposed to thawing epi- 
sodes during sample transport from the remote location, 
thus we predicted a skewed microbial community compo- 
sition in comparison with EtOH-stored gut samples and 
tested the extent and constancy of this effect. 


Methods 
Study site and sampling methods 


Tadpoles of Nanorana parkeri in Gosner stages (GOSNER 
1960) 26-30 were collected on 2 July 2018 from the cen- 
tral Qinghai-Tibetan Plateau, China, in a small pond di- 
rectly nearby a small tributary stream, in the area of Lake 
Nam Co (30.82840° N, 91.06397° E; 4,740 m above sea lev- 
el). Tadpoles were collected using dip nets and placed indi- 
vidually into sterile Whirl-Pak bags together with a small 
amount of clean water from the respective water bodies. A 
pond water control sample was collected by dipping a ster- 
ile swab into the water, placing it into a cryotube and freez- 
ing it (list of taxa found from pond water are outlined in 
Supplementary Table S1). Specimens were kept overnight 
(ca. 14 h) in the bags, anesthetized with tricaine methane- 
sulfonate (MS222; Sigma-Aldrich, St. Louis, MO, USA) so- 
lution and subsequently overdosed using MS222. Sacrific- 


ing of the specimens was reported to the animal welfare 
officer of Technische Universitat Braunschweig in compli- 
ance with German animal welfare law, section 4(3). We col- 
lected the feces accumulated in the bags using pipettes, as 
well as a portion of hindgut of every individual. In sum- 
mary, we obtained three replicate samples per specimen: 
(1) feces, stored in 96% EtOH; (2) hindgut, stored in 96% 
EtOH; (3) hindgut, frozen at -20°C right after collection 
(Fig. 1). Samples were frozen upon collection, but as typical 
for suboptimal fieldwork conditions, underwent two thaw- 
ing-freezing cycles during transport to the lab, with tem- 
peratures of 10-15°C for a limited time (< 12 h each time). 
In the laboratory, all samples were stored at -20°C until 
further processing. 


Molecular analysis 


Following the manufacturer’s instructions, DNA was ex- 
tracted using a DNeasy PowerSoil Kit (QIAGEN, Germa- 
ny). PCR was performed using the forward primer 515F 
(5 -GTGCCAGCMGCCGCGGTAA-3’) and reverse prim- 
er 806R (5 -GGACTACHVGGGTWTCTAAT-3’) to target 
the V4 region of 16S rRNA gene (Caporaso et al. 2011). 
The used primer and tag (unique molecular identifier for 
multiplexing) combinations for each sample are specified 
in Supplementary Table S2. For PCR, the 25 ul mixture per 
sample comprised of 2 ul DNA (3 ul for repeated samples), 
0.5 ul each of the primer (10 uM), 4 ul 5x HOT FirePol® 
Blend Master Mix (Solis BioDyne, Tartu, Estonia) and the 
rest of the volume was filled with sterile dH,O. PCR was 
carried out in two replications in the following thermocy- 
cling conditions: an initial 15 min at 95°C, followed by 35 
cycles of 94°C for 45 s, 50°C for 60 s, 72°C for 90 s, and a fi- 
nal cycle of 10 min at 72°C. PCR products per sample were 
pooled and their relative quantity was estimated during 
gel electrophoresis of 5 ul DNA sample on 1% agarose gel. 
Based on gel band intensity, all PCR products were pooled 
at approximately equimolar concentration. The DNA li- 
brary was purified using Favor-Prep™ Gel/PCR Purifica- 
tion Kit (Favorgen-Biotech Corp., Vienna, Austria), fol- 
lowing the manufacturer’s instructions. Sequencing was 
performed on an Illumina MiSeq instrument (2x250) us- 
ing the MiSeq Reagent Kit v2. Steps of DNA extraction, 
PCR and sequencing included both negative and positive 
controls (Supplementary Table S2). Sequencing data (in- 
cluding control samples) have been deposited in the Se- 
quence Read Archive (SRA): BioProject ID PRJNA533915. 


Bioinformatics 


The paired-end sequence data was processed in QUME 
(v1.9.0; CAPORASO et al. 2010) using the Environmen- 
tal Microbiome and Bioinformatic Analysis Platform of 
School of Public Health in Lanzhou University. Data anal- 
ysis methods were as described previously (Lı et al. 2016b). 
Briefly, paired-end sequences were joined using Flash soft- 
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ware (v1.2.8; MAGOC & SALZBERG 2011). Those sequences 
with length < 300 bp, average base quality score < 30 or 
ambiguous bases, were removed for the downstream anal- 
ysis. The Uchime algorithm (EDGAR et al. 2011) was used to 
filter out potential chimeric reads. The filtered sequences 
were clustered into operational taxonomic units (OTUs) 
at a 97% identity threshold using UCLUST algorithm 
(EDGAR 2010). Taxonomy was assigned using the Ribos- 
omal Database Project classifier (WANG et al. 2007). OTUs 
not classifying to Bacteria, Eukaryota and Archaea line- 
ages, were removed. Although Archaea may represent an 
important part of gut microbiome (RAYMANN et al. 2017), 
we removed this group from the analyses because of the 
low specificity of the used primers (515F/806R) to Archaea 
(BAHRAM et al. 2018a, RAYMANN et al. 2017) (our data set 
contained < 0.001% of Archaea sequences). To account for 
the unequal sequencing depth, each sample was rarefied to 
the same number of reads (5,654 sequences). The latter led 
to discarding one of the frozen gut sample from the down- 
stream analysis. The OTU table was further filtered to re- 
move singleton OTUs and low-abundance read records of 
OTUs per sample (< 10 reads). After these steps, cyanobac- 
terial OTUs were removed from the analysis as these taxa 
likely do not represent true gut microbiota (e.g. Nostoco- 
phycideae, Synechococcophycideae, Oscillatoriophycideae 
(PATZELT et al. 2014)). Additionally, taxa that were assigned 
to ‘mitochondria, were also excluded from the analyses as 
potential mitochondrial reads of the host itself. The filtered 
OTU table used for the analyses is specified in Supplemen- 
tary Table S3. 


Statistical analysis 


The effect of sample type (EtOH-stored feces, EtOH-stored 
gut and frozen gut) on log-transformed OTU richness and 
phylogenetic diversity (PD) were tested using Kruskal- 
Wallis analysis. OTU richness correlation analyses be- 
tween samples types were performed using Spearman rank 
correlation analysis. The effect of these factors on the bac- 
terial OTU composition was analyzed using PERMANO- 
VA+ (ANDERSON 2005) with 9,999 permutations (Type II 
SS) in PRIMER v6 (CLARKE & GORLEY 2006). For the 
PERMANOVA analysis and non-metric multidimension- 


Feces 
Gut 


al scaling (NMDS) graphs, we used Hellinger-transformed 
Bray-Curtis as well as UniFrac distance (unweighted) OTU 
matrices. UniFrac distance and PD were calculated by ap- 
plying the PhyloMeasures package (v2.1; TstROGIANNIS 
& SANDEL 2016) in R (R Core-Team 2020) using a Maxi- 
mum-Likelihood based phylogenetic 16Sr RNA gene tree 
generated with RAxML (Stamataxis 2014). The similar- 
ity between Hellinger-transformed Bray-Curtis distance 
based NMDS ordinations of each sample type was test- 
ed using Procrustes analyses using the R-package ‘vegan’ 
(OKSANEN et al. 2015). We used indicator species analysis 
(‘indicspecies’ v1.7.6 package in R; DE CACERES et al. 2016) 
to determine which OTUs are significantly (using 9,999 
permutations) associated with particular sample types. A 
Venn diagram of shared OTUs between sample types was 
created using mothur (ScHLOss et al. 2009). 


Results 


Our analysis comprised samples of 12 individual tadpoles 
(Fig. 1). Replicate samples of EtOH-stored feces, EtOH- 
stored gut and frozen gut per specimen demonstrated sig- 
nificant (Kruskal-Wallis test: P < 0.001) and high (3.8-fold) 
OTU richness variability (Supplementary Fig. S1). This 
pattern was driven by the frozen gut samples that exhibited 
significantly lower OTU richness and phylogenetic diver- 
sity (PD) values as compared to the other samples (Figs 2a- 
b). Consequently, there was no strong positive sample by 
sample correlation of the detected OTU richness between 
frozen gut and other treatments (Spearman R = 0.179, P = 
0.599 for gut (EtOH) vs. gut (frozen); Spearman R = -0.378, 
P = 0.251 for feces (EtOH) vs. gut (frozen); Supplementary 
Fig. S2). This correlation was highest between EtOH pre- 
served samples (gut vs. feces; Spearman R = 0.494), how- 
ever, the association was statistically non-significant (P = 
0.103; Supplementary Fig. S2). 

The communities of detected microbiota were signifi- 
cantly different between sample types (PERMANOVA 
P < 0.001; Table 1; Figs 2c-d; Fig. 3). EtOH-stored sam- 
ples vs. frozen gut samples demonstrated to harbor low 
proportions of shared OTUs, 14.1% and 15.6% (Fig. 4) for 
EtOH-stored gut vs. frozen and EtOH-stored feces vs. fro- 
zen gut samples, respectively. Accordingly, based on the 


EtOH (96%) 


room temp. 


SAMPLE TRANSPORT 


ANALYSES 


Figure 1. Illustration of the study design. Three treatments: ethanol (EtOH) stored feces and gut, and frozen gut samples of Nanorana 
parkeri. Note that the frozen samples went through thawing episodes during sample transport from the field to the laboratory. 
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Table 1. PERMANOVA results for Nanorana parkeri samples 
with factor sample type (EtOH-stored feces, EtOH-stored gut 
and frozen gut). PW denote pairwise comparison between EtOH- 
stored gut and frozen gut (PW1), EtOH-stored feces and fro- 
zen gut (PW2), EtOH-stored feces and EtOH-stored gut (PW3) 
samples. 


OTU matrix R? ä P P PWI P PW2 P PW3 
Bray-Curtis 0.470 <0.001 <0.001 <0.001 <0.001 
Unifrac 0.361 <0.001 <0.001 <0.001 0.007 


Procrustes analyses, the pairwise comparisons of the lat- 
ter sample types showed the highest differences (i.e. weak- 
est Procrustes correlation) with the correlation values of 


0.346 (P = 0.474) and 0.370 (P = 0.370) for EtOH-stored 
gut vs. frozen gut and EtOH-stored feces vs. frozen gut 
samples, respectively. The highest similarity, but margin- 
ally non-significant Procrustes correlation was observed 
between EtOH-stored feces and EtOH-stored gut samples 
(Procrustes correlation = 0.523, P = 0.072). The proportion 
of shared OTUs between the latter sample types was 26.2% 
(Fig. 4). 

Indicator species analysis revealed a large number of 
OTUs that were differentially abundant in across samples; 
13 OTUs for EtOH-stored feces, 28 OTUs for EtOH-stored 
gut (30 OTUs for the combination of EtOH-stored sam- 
ples) and 12 OTUs for the frozen gut samples (Supplemen- 
tary Table S4). The combination of EtOH-stored samples 
(feces + gut) harbored many indicator OTUs, whereas 
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Figure 2. a) OTU richness (log transformed) and b) phylogenetic diversity (log transformed) for EtOH-stored feces, EtOH-stored 
gut and frozen gut samples of Nanorana parkeri. Whiskers denote 95% confidence intervals. Different letter combinations on top of 
whiskers denote significant difference between groups (a < 0.05) as based on Tukey HSD test. c, d) NMDS graphs as based on Bray- 
Curtis similarity and UniFrac distance of bacterial community composition. e) Bar plots showing relative abundance of reads assigned 


to bacterial phyla, summarized for each treatment. 
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no OTUs were assigned to be specific for sample combi- 
nations with frozen gut samples (i.e. EtOH-stored feces + 
frozen gut, EtOH-stored gut + frozen gut) (Supplementary 
Table S4). This further suggests an overall higher similar- 
ity of the two types of EtOH-stored samples, which is also 
evident from the NMDS graphs (Figs 2c-d). 

Frozen gut samples showed a consistent increase of 
Proteobacteria (sequence abundance; Fig. 2e), which was 
mainly caused by ten of the indicator OTUs, all belong- 
ing to the Gammaproteobacteria: two undetermined 
Aeromonadaceae, two undetermined Enterobacteriace- 
ae, two Klebsiella (Enterobacteriaceae), two Pseudomonas 
(Pseudomonadaceae), and two Shewanella (Shewanel- 
laceae). These OTUs represented a very high proportion 
of reads in the frozen gut samples (90.1%), but most of 
them were completely absent from all of the other sample 
types. EtOH-stored gut samples, however, showed high- 
est relative abundance of Firmicutes (sequence abundance; 
Fig. 2e). This pattern was particularly driven by Clostridia 
of which 25 OTUs were identified as indicators, and thus 
were relatively more abundant in the gut (EtOH) samples 
(Supplementary Table S4), forming 26.4% of the reads on 
EtOH-stored gut samples. 


Feces (EtOH) Gut (EtOH) 


Relative abundance of OTUs 


Specimen 1 
Specimen 2 
Specimen 3 
Specimen 4 
Specimen 5 
Specimen 6 
Specimen 7 
Specimen 8 
Specimen 9 
Specimen 10 
Specimen 11 
Specimen12 
Specimen 1 
Specimen 2 
Specimen 3 
Specimen 4 
Specimen 5 
Specimen 6 
Specimen 7 
Specimen 8 


Figure 3. Relative abundance of bacterial OTUs in all samples. 


100 


1.00% 
0.75. 
0.50- 
0.25. 
0.00- 


Specimen 9 


Discussion 


As the number of DNA sequences of environmental or 
host-associated samples in public databases increase and 
acquiring such data becomes a routine approach in micro- 
bial ecology, meta-analyses of comprehensive “big data” 
sets is becoming a promising research direction, leading 
to important insights into general patterns of bacterial di- 
versity (ADAMS et al. 2015, GONZALEZ et al. 2018, ROCCA et 
al. 2018, KUENEMAN et al. 2019). However, considering the 
many technical factors influencing the outcome of ampli- 
con analyses, including sample preservation, DNA extrac- 
tion, PCR conditions and sequencing methods (GOODRICH 
et al. 2014, CLOONEY et al. 2016, ANSLAN et al. 2018), it is of 
high importance to ascertain that such meta-analyses in- 
deed recover biological patterns and not methodological 
differences among studies. 

Our study exemplifies that the recovery of host-asso- 
ciated microbiota richness and community structure may 
vary among sampling substrates. The two EtOH-stored 
substrates, gut content recovered by dissection and fecal 
samples, overall, revealed rather consistent patterns of 
OTU richness and phylogenetic diversity. However, con- 
sidering community composition, important differenc- 
es between gut and fecal samples were apparent, where 
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only about one fourth of the OTUs was shared between 
these sample types (Fig. 4). Similarly, ZHov et al. (2020) 
recently reported that the microbial composition from fe- 
cal samples of adult cane toads (Rhinella marina) were sig- 
nificantly different from the intestine (gut) samples. Simi- 
lar results were found by GRIFFIN et al. (2020) by studing 
blue mussels (Mytilus edulis), where they discourage the 
use of fecal samples to explore the gut microbial commu- 
nities. One factor affecting the community difference be- 
tween gut and feces could be the exposure to oxygen upon 
leaving the intestinal tract (RAGo et al. 2017). Although 
we included a pond water control in our study to account 
for external contamination, the exposure to environmen- 
tal contamination of fecal samples may further have fa- 
cilitated the changes in microbial composition associated 
with feces samples (SONG et al. 2018, GRIFFIN et al. 2020). 
In this study and previous studies by ZHov et al. (2020) 
and GRIFFIN et al. (2020), a single DNA extraction meth- 
od was applied, but by utilizing two different DNA extra- 
cion methods to study the gut microbiota of the Asiatic 
toad tadpoles (Bufo gargarizans), SONG et al. (2018) report- 
ed contrasting results. The results from samples subjected 
to TIANamp stool DNA kit demonstrated significant dif- 
ferences between gut and feces, but interestingly, samples 
subjected to phenol-chloroform DNA extraction exhibited 
consistent microbial community composition patterns be- 
tween these two sample types. Based on these results, the 
authors suggested that feces of amphibian tadpoles may be 
effectively used to study their gut microbiota. This high- 


Gut (EtOH) 


156 
44.8% 


Figure 4. Venn diagram of shared OTUs between sample types. 
Total number of OTUs is 348. The number of OTUs in feces 
(EtOH), gut (EtOH) and gut (frozen) samples is 158, 253 and 88, 
respectively. The number of shared OTUs between feces (EtOH) 
and gut (EtOH) is 91 (26.2%). The number of shared OTUs be- 
tween feces (EtOH) and gut (frozen) samples is 54 (15.6%), and 
49 (14.1%) between gut (EtOH) and gut (frozen) samples. Total 
numbers of shared OTUs is 43 (12.4%). 


lights also the effect of DNA extraction method on the re- 
covered microbial composition, which has been reported 
in several other studies (e.g. CHAPAGAIN et al. 2019, FIE- 
DOROVA et al. 2019). Althogether, it is obvious that merg- 
ing results based on these different types of substrates (gut 
and feces) for a meta-analyses may introduce biases, and 
thus they should be combined with great caution, and only 
when large differences between hosts are expected (which 
then would override the influences of methodology and 
substrate). 

An even more divergent pattern was found between 
commonly used sample preservation methods (EtOH 
vs. freezing), with enormous differences both in bacteri- 
al richness and community structure. Consistently, in al- 
most all individual tadpole (frozen) gut samples, Proteo- 
bacteria had enormously increased relative abundances, 
whereas relative abundances of Firmicutes and Verruco- 
microbia were much lower (Fig. 2e; Fig. 3) compared to the 
EtOH-preserved samples. As summarized by Kout (2017), 
blooming of certain taxa can change the composition of 
gut or fecal bacterial communities (CHOO et al. 2015, BECK- 
ERS et al. 2017). We hypothesize this is what happened in 
our samples upon thawing periods during sample trans- 
port. The blooming hypothesis is also supported by the fact 
that the Proteobacteria increase was caused by a limited 
number of bacterial OTUs, and most strongly influenced 
by only seven OTUs. Moreover, these OTUs were totally 
absent from the EtOH-stored Nanorana gut samples. How- 
ever, several studies on soil, human- and insect-associated 
microbiomes have revealed that differences among preser- 
vation methods are smaller than those between taxa and 
individuals, thus validating meta-analyses based on differ- 
ently stored samples (LAUBER et al. 2010, DOMINIANNI et 
al. 2014, HAMMER et al. 2015, BLEKHMAN et al. 2016, SONG 
et al. 2016). A study about sample preservation methods of 
fecal microbiota of spider monkeys revealed that the mi- 
crobial community composition of EtOH-stored and fro- 
zen samples were similar to fresh ones (HALE et al. 2015). 
Thus, it is expected that the latter sample storing methods 
are producing comparable results. Furthermore, a study on 
insect-associated microbiomes has suggested that the sam- 
ple storage method (freezing, ethanol, dimethyl sulfoxide, 
cetrimonium bromide, storage without any preservative) 
has no or minor effect on microbiome composition (HAM- 
MER et al. 2015). However, under typical field sampling 
conditions, such as in this study, continuous deep-freez- 
ing of samples cannot be always ensured, which implies the 
possibility of radical effects on gut microbiomes. Sampling 
ungulate feces in the wild, MENKE et al. (2015) observed 
only moderate shifts of the microbiome during 2-4 days 
but radical changes afterwards, usually following rain. On 
the contrary, BECKERS et al. (2017) observed an important 
decrease in bacterial diversity in horse feces already after 
approximately 4 hours. A significant decrease of Bacteroi- 
detes has been reported from fecal samples (of humans and 
monkeys) exposed to room temperature or to natural en- 
vironmental conditions (ROESCH et al. 2009, HALE et al. 
2016), which we also found to be the case for the frozen 
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samples that were exposed to thawing in current study. 
Our study confirms that in certain cases, a rise of “bloom” 
bacteria can completely obscure the original microbiome 
composition in the samples, and several bacterial families 
might be particularly prone to contain such rapid growth 
taxa; for instance, Enterobacteriaceae and Pseudomonada- 
ceae were associated with microbiome shifts both in this 
study and in that of BECKERS et al. (2017). 


Conclusion 


Our case study confirms that both substrate (gut content 
vs. feces) and preservation method can influence the anal- 
ysis of intestinal microbiomes, and provides an example 
from tadpoles. Differences between substrates and preser- 
vation methods are here shown for samples from exactly 
the same individuals, sampled at the same time point, thus 
excluding these factors that might influence microbiome 
structure. Although non-invasive collection of feces may 
be preferred over dissection for ethical reasons, it is im- 
portant to acknowledge the recovery of potentially differ- 
ent microbiome communities between these treatments, 
as demonstrated in this study. When using fecal samples, 
one must be also aware of environmental contaminations 
of these samples, which could be minimized by collecting 
control samples in the immediate environment. Although 
the latter was attempted in our study, significant differ- 
ences between microbial composition of gut and feces 
were still obvious. Emphasizing the thawing cycles of fro- 
zen samples, the strongest differences were found between 
preservation methods and demonstrate that blooming of 
contaminant taxa can completely distort the bacterial com- 
munity in samples of intestinal microbiome of aquatic ver- 
tebrates, within only a few hours of thawing as it is com- 
mon under field conditions. 
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